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^ : ABSTRACT 

(N , 

[V] ' We carry out a spectral analysis of the archival FUSE spectrum of the VY 

^ ■ Scl nova-like cataclysmic variable MV Lyrae obtained in the high state. We find 

"-^ ' that standard disk models fail to fit the flux in the shorter wavelengths of FUSE 

I ' (^ < 950A). An improved fit is obtained by including a modeling of the boundary 

^-^ • layer at the inner edge of the disk. The result of the modeling shows that in 

52 ■ the high state the disk has a moderate accretion rate M ?« 2 x lO~^M0/yr, 

'""'■ a low inclination, a boundary layer with a temperature of ~ 100,000K, and 

T-H ■ size ~ Q.2QRwd, and the white dwarf is possibly heated up to a temperature 

^'- T ^ 50, OOOK or higher. 

a^: 

t:;^- ■ Subject headings: accretion, accretion disks - novae, cataclysmic variables - white 

[--. ■ dwarfs 

o; 

^H , 

^ ■ 1. Introduction 

X ■ 

^ : 1.1. The VY Scl Nova-like System MV Lyrae 

Cataclysmic variables (CVs) are evolved compact binaries in which the primary, a white 
dwarf (WD), accretes matter and angular momentum from the secondary, a star approxi- 
mately on the main sequence filling its Roche-lobe. In dwarf nova (DN) and non-magnetic 
nova-like (NL) systems (two classes of CVs), the matter is transferred by means of an accre- 
tion disk around the WD. DN systems are found mostly in a faint, quiescent state with a low 
mass transfer rate (M ~ lO^"*^^ — IQ^^^Mq/ji), and every few weeks to months they transit 
into a bright, outburst state lasting days to weeks, as the mass transfer suddenly increases 
by several orders of magnitude (M ~ 10"^ — lO'^M^^/yr). The transition to outburst in DNs 



is believed to be due to a disk instability ( IShafter et al. Ill986l : ICannizzo Ill988l . Il998l ). NL 
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systems are found mainly in a high state phase, characterized by a large mass accretion rate, 
and from time to time the mass accretion rate suddenly drops by several orders of magni- 
tude. The transition from high state to low state in NL is caused by a throttling of the mass 
transfer, possibly due to magnetic activity of the secondary star (see IWarner I (Il995l ) for a 
review of CVs). With time, the accreted matter in CV systems forms a hydrogen-rich layer 
that covers the surface of the WD. When enough material is accreted (every few thousand 
years) , the hydrogen-rich material undergoes a thermonuclear runaway and the system forms 
a classical nova event. 

MV Lyr is a VY Scl nova-like system, spending most of its time in the high state 
and occasionally undergoing short- duration drops in brightness. When this happe ns, the 
magnitude of MV Lyr drops from ^ si 12 — 13 to \^ ~ 16 — 18 fJHoard et al. 1 120041 ). From 



the AAVSO online data (http://www.aavso.org/) it appears that MV Lyr can stay in the 
high state for up to ~ Syr, and then for a period of a few years to ~ lOyr it spends its time 
alternating between high state and low state on a time scale of a few months to a year. The 
orbital period of the b inary in MV Lyr is P = 3.19hr and the mas s ratio was found to be 
q = M2nd/M^d = 0.4 (ISchneider et al. Ill98ll : ISkillman et al. Ill995l ). Because of the small 



radial velocity amplitudes and the lack of eclipses, the inclination is con strained to be in the 
range i = 10° — 13°, though it has been remarked ( iLinnell et al. Il2005[ ) that the inclination 
could be as low as z = 7°. 

Observations of MV Lyr during a low state with t he International Ultraviol et Explorer 
(WE) indicated a hot W D possibly reaching 50,000K (ISzkody &: Downes Ill982[ ) or higher 
( IChiappetti et al. Ill982[ ). These observations were later confirmed when MV Lyr was ob- 
served with FUSE on July 7, 2002, during a low state that lasted for about 8 months. 
Hoard et al. I ( 120041 ) carried out a spectral analysis of the FUSE spectrum and found that 
the WD must have a temperature of 47,000K, a gravity logg = 8.25, a projected rotational 
velocity Vrot sin i = 200km/s, subsolar abundances Z = 0.3 x Zq, and a distance of 505±50pc. 
At the time of the observations, the magnitude of MV Lyr was po ssibly as low as V ~ 18 



with a mass accret ion rate no niore th an M ^ 3 x 10 ^^Mq /yi (JHoard et al. Il2004f ). A 



follow up analysis ( ILinnell et al. Il2005[ ) was carried out to study the different states of MV 
Lyr and its secondary, based on archival lUE spectra, a HST/ffr/ ^* snapshot and optical 



spectra. Using the parameters they derived in iHoard et al. I ( 120041 ). ILinnell et al. I (|2005[ ) 



found that standard disk models did not provide adequate fits to the spectra; instead, the 
disk had to be modified by, for example, a truncation of the inner disk, an isothermal radial 
temperature profile in the outer disk, or both. These modified disk models provided some 
satisfactory fit to the spectra obtained in the high state, leading to a mass accretion rate of 
the order of 3 x lO^^Mo/yr. 
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In the present analysis, the standard disk model fails again to provide a good fit to the 
observed {FUSE) spectrum of the system in the high state. This, however, points to a more 
general problem, namely that the standard di sk model has failed to fit the observed spectra 
of many other CVs observed in the high state (jPuebla et al. 1120071 ) . The direct implication is 
that the standard disk model itself is inadequate in the ultraviolet (UV). It is in this context 
that we analyze here the FUSE spectrum of MV Lyr in its high state. For this reason, we 
review next the problems encountered in the modeling of CVs in the high state. 



1.2. UV Spectral Modeling of CVs in the High State 

In high state systems, the UV emission predominantly originates from the disk, and 
for that reason these systems became of special interest, as they were expected to have an 
accretion disk radial temperature profile given by the arialytica l expression Teff{r) of the 
standard disk model. More than two decades ago, IWade I (Il988l ) used lUE spectra to show 
that CV NL systems (high state systems) systematically disagree with the standard disk 
model when either blackbody spectra or Kurucz stellar model spectra were use d to re pre- 
sent the accretion disk. A much improved version was developed by iHubeny I ( 119901 ) for 
modeling annuli of accretion disks that explicitly included calculation of synthetic spectra 
using the standard disk model (the c odes TLDISK/TLUSTY , SYNSPEC and DISKSYN), 



and sp ectra of disks were computed (IWade &: Hubenv 



19981) a s suming the st andard disk 



model (jShakura fc Sunyaev Ill973l : iLynden-Bell fc Pringle Ill974l : iPringle Ill98ll ). UV spec- 
tra (HST, HUT, lUE, FUSE) of DNs in outburst and high brightness states of NLs were 
successfully modeled for a number of systems with synthetic spectra of optically thick disks 
(iKnigge fc Long Il2002[ [Hamilton et al. \\20Ql\ \. 



However, with time, it became cle ar that a significant number of systems still could not 



be fit ted using the standa rd disk rnodel (ILong et al. lll994 : lKnigge et al. 1119971 : iRobinson et al. 



19991 ) as first shown by IWade I ( 119881 ). lOrosz fc Wade 



( 120031 ) suggested a disk tempera- 



ture profile {T{r)) fiatter than the standard disk model, to remediate at least partially to 
the problem, invoking a possible non-steady mass accretion rate and/or irradiation in the 
outer part of the disk. More rec ently, the UV spectra of RW Sex, V3885 Sgr, & UX UMa 
( jLinnell et al. I l2008l . l2009l . l2010l ) were also modeled using a modified temperature profile, 
and it was suggested that the outer disk temperature might increase due to its interaction 
wi th the impacting LI stream. Maybe the most significant of all is the statistical study 



of iPuebla et al. I ( 120071 ) of 33 disk-systems which indicates that a revision of the standard 



model temperature profile, at least in the innermost part of the disk, is required. 



These studies concentrated mainly on using UV spectra (e.g. lUE, RST / STIS, FUSE), 
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and the contribution to the UV continuum conies mainly from regions in the system where 
the temperature ranges roughly between ~10,000K and ~100,000K. Lower temperatures 
peak in the optical, while higher temperatures peak in the extreme UV. The WD star and 
the disk are the main contributors to the UV flux, while the secondary star and the bright 
spot (the edge of the disk impacted by the incoming stream of matter) peak in the optical 
and near UV, respectively. For disks in CVs with mass accretion rate ~ lO~^M0/yr, the 
inner disk reaches a temperature higher than 50,000K. These disk models, therefore, have a 
significant UV flux contribution from the inner region of the disk. 

That same innermost part of the disk is the region where the matter in the disk has to 
decelerate from its Keplerian motion in order to be accreted onto the more slowly rotating 
stellar surface, as the centrifugal force keeps the matter from accreting. This region is known 
as the boundary layer (BL) and up to half the accretion luminosity can be released in this 
final phase of the accretion process. Therefore, the emission and heating of the BL have to be 
taken into account when mode l ing an accretio n disk around a non- (or weakly- ) magnetic star 
(JLynden-Bell fc Pringle Ill974j : iPringle 1119811 ). It has now become clear that the temperature 
profile in the inner disk is affected by the heating from the BL. Most importantly, BLs in 



high state have a temperature of the order o f 100,000K-200,000K (jPopham fc Narayan Ill995 



Go don et al. I Il995l : iPiro fc Bildsten 1120041). which is much lowe r than the first analytical 
estimates of ~500,000K (IPringle 1119771 : IPringle &: Savonije Ill979l ). Consequently, the BL of 
CV disk systems in the high state contributes significantly to the UV, even though it was 
predicted to emit only in the soft X-ray range (see also next section). Therefore, the BL 
has to be included in the disk modeling and the inner disk temperature profile has to be 
modified accordingly. 

It is the purpose of the present work to include an elementary (and preliminary) model 
of the BL in the spectral modeling of the FUSE spectrum of MV Lyr taken during its high 
state. Since the BL has an elevated temperature, its contribution will be more pronounced 
in the shorter wavelengths of the spectral range of FUSE. 

In the next section a brief review of the standard disk and BL models are given, and it 
is shown that the theory does indeed predict a BL contributing a non- negligible flux in the 
far UV. The FUSE spectrum of MV Lyr together with the details of the synthetic spectral 
modeling are presented in section 3. The results follow in section 4, and are discussed in the 
concluding section. 
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2. The Star-Disk Boundary Layer 



2.1. The Standard Disk Model 



The total luminosity available through accretion is given by: 

_GMj£ 

-'-'ace T-, ) 

-fX* 



:ii 



where G is the gravitational constant, M is the mass accretion rate, M^, is the mass of the 
accreting star and i?* is its radius. This expression simply reflects the release of gravitational 
potential energy of the material brought onto the surface of the star {r = R^) from inflnity 
(r — )■ cxd) at a rate M (mass/time). 



In the standard disk theory (jShakura fc Sunyaev Ill973l : iLynden-Bell fc Pringle Ill974j ). 
the accretion disk is axisymmetric, and geometrically thin in the vertical dimension: it has 
a vertical thickness H such that H/r << 1. The disk is treated with a one-dimensional 
radial approximation, in which the radial pressure is negligible and the matter in the disk 
is assumed to be rotating at Keplerian speed. The energy dissipated by viscous processes 
between adjacent rings of matter is instantly radiated locally in the vertical z-direction. The 
disk total luminosity Ldisk is half the accretion luminosity 



-'disk 



GM^M 

2Ri: 



(2) 



and each face of the disk radiates Ldisk/'2. This is so, because the matter at the inner edge 
of the disk {r k. R^) retains the remaining accretion energy in the form of kinetic energy, as 
it rotates at Keplerian speed. 



The effective surface ternperature of suc h a disk is given by (IShakura fc Sunyaev Ill973 
Lvnden-BeU k Pringle Ill974 



Pringle I Il981[ l: 



eff[ 



X 



1/2x1/4 



where x = r/R^ 
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and a the Stefan-Boltzmann constant. For practical purpose the last relation is usually 

written 
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From this it is seen that accretion disks around WD stars emit their energy in the optical and 
UV. The above relation is obtair ied by assumin g; a no-shear boundary condition, dVt/dr = 0, 
at the stellar surface r = R^ ( iPringle I Il98ll ). Therefore, this model does not take into 
account the relatively slow rotation of the WD, and it gives a disk temperature T=OK (!) 
at the stellar surface (r = R^), and a maximum temperature Tmax = 0.488To at x = 1.36. 
This model also neglects the spin up of the star by the disk. 

This relation, however, gives a good approximation of the effective surface temperature 
in the disk at larger radii (r >> R^). For a typical high mass accretion rate of M = 
1 xlO~^M0/yr, the maximum temperature in the inner disk reaches about 50,000K-100,000K, 
and the inner disk become the strongest emission source of the system in the UV. During the 
low state of NLs (or quiescent state of DNs) the mass mass accretion rate decreases down to 
M ^ 10^^'^Mq/jy, while the maximum temperature drops to less than f«10,000K, making 
the disk peak in the optical and barely emit any flux in the UV. The standard disk model 
is a good approximation for r >> i?*, but in the inner disk region, one needs to model the 
BL between the slowly rotating stellar surface and its Keplerian accretion disk. We would 
like also to note here that contrary to accreting neutron stars, the irradiation of the disk by 



the heated WD and BL is negligible in accreting WDs (Ivan Paradijs fc McClintock I Il994 



Shahbaz fc Kuulkers Ill998l : King lllQQSl ). and one does not need to take it into account when 
modeling accretion disks in CVs. 



2.2. The Boundary Layer 

A ccreti ng WDs in CVs have stellar rotational velocities of the order of a few lOOkm/s 

(ISion 1119981 ) or about 10% of the Keplerian velocity at one stellar radius (a few lOOOkm/s). 
Consequently, the theory (JLynden-Bell &: Pringle I Il974j : iPringle I Il98ll ) predicts that the 
material at the inner edge of the disk rotating at Keplerian speed has to adjust itself to 
the slowly rotating stellar surface and, therefore, dissipates its remaining rotational kinetic 
energy. The remaining rotational kinetic energy liberated in the BL {Lbl) is nearly equal 
to half of the to tal gravitational energy of the accreting matter {Laccf^ = Ldisk), namely 
kluzniaklll987l l: 

Lbl = Ldisk ( 1 - ^ /^ x ) • (6) 



Because of its small radial extent and the large amount of energy dissipated there, the 
BL was expected to emit only in the X-ray bands. At high mass accretion rates (M ~ 
10~^ — 10~^ MQyr~^) the BL should be o ptically thick and its temperature was pre dicted to 
be in the range 200,000K to 500,000/s: (jPringle I [l977l : IPringle fc Savoniie I[l979[ l. At low 
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mass accretion rates, the BL should be optically thin and emit in the hard X-ray (20keV), 
due to either strong shocks or the formation of a X-ray emitting corona around the WD 
dPrinde fc Savoniie lEoToi king fc Shaviv Ill984[ ). 



The first one-dimensional cornputations of the BL ( Regev I Il983l : iPapaloizou fc Stanley 



19861 : [Stanley fc Papaloizou I Il987l: iRegev fc Hougerat 



19881 ) were carried out by assum- 



ing the "shm disk equations" (ISzuszkiewicz et al. Ill988l ): the equations are written and 
solved for vertically averaged quantities but also include the transport of energy (radiation 
and advection) in the radial dimension. The results from these one- dimensional computa- 
tions were in general agreement with the predictions and showed that the BL should be 
narro w with a temperature increasing rapidly inward. However, rnore realistic ID simul a- 
tions ( Popham fc Naravan Ill995l : ICodon et aLlll995l : ICollins Ill997[ ICoUins et aLlll998al Jb. 
2000a| ]bl) showed that the BL extends radially and vertically more than expected. The result 
was that a t high accretion rate s the BL is not as hot as expected, with a temperature around 
125,000K (ICodon et al~lll995h to 300,000K JPopham fc Naravan |[l995l ). 



Numerical works in 2D ( iRobertson fc Frank Ill986l : iKley fc Hensler II 19871 : iKley lll989al Jb 



199ll ) did not clearl y reveal w 
e.g. as claimed by iDurisen I ( 



over the entire WD surface (JMacDonald 



lethe r the accreting material accuni u lates in an equatorial belt 
1977) and iKippenhahn &: Thomas Id 19 78[) ') or spreads quickly 



1983 



using the semi-analytical shallow- water approximation de veloped by llnogamov fc Sunyaev 



Livio fc Truran Ill987l). But more recently. 



( 119991 ) for accreting neutron stars, iPiro &: Bildsten I ( 120041 ) found that at high mass accretion 
rate (M > lO^^g s~^ ~ 1.6 x 10~^MQ/yr) the combined effect of centrifugal force and 
pressure gives rise to two latit udinal rings of enhanced brigh tness above and below the WD's 
equator: the "spread layer" (llnogamov fc Sunyaev I Il999l ). At lower mass accretion rate 
(M < lO^^g s^^ ^ 1.6 X lO~*M0/yr) the spreading is negligible and most of the dissipated 
energy is radiated back into the accretion disk. The effective temperature (2 x lO^K-5 x lO^K) 
and rotational velocity (1 x lO^cm/s-3 x 10*cm/s) of the spread BL are similar to those of 
the one- dimensional BL models (at least where they overlap in the para meter space). Some 



more recent simulations ha ve been carried out to study the spread layer (JFisker fc Balsara 



20051 : iBalsara et al. II2009I ). however, these two-dimensional simulations do not yet include 
radiation hydrodynamics. 

The significance of a cooler BL means that its emission in the UV band is stronger 
and cannot be neglected when computing the luminosity of the disk and WD. In the next 
subsection (sec. 2. 3) we summarize the BL emission one may expect to observe in the UV at 
different mass accretion rates based on the above mentioned investigations. 



2.3. The Emission of the Boundary Layer 



For M Ri 2 X lO~^M0/yr and higher, one expects the BL to spread onto the stellar 
surface, with a temperature of the order of 2 x lO^K-5 x lO^K, with a velocity of the order 
of 0.3 — 1 X Vf^(R^). The spread layer should be observable above and below the disk 
(jPiro fc Bildsten Il2004f ). At such temperatures the BL emits mainly in the soft X-ray with 
some emission in the extreme UV, and possibly also far UV. The remaining part of the WD 
not covered by the spread layer emits in the UV and is likely to have an elevated temperature. 
The disk itself emits in the UV (far UV to near UV as the radius increases) and optical 
(outer disk). As the mass accretion rate increases, the thickness of the boundary/spread 
layer increases and the matter becomes advection-dominated as the heating rate surpasses 
the cooling rate. In this extre me case, whi ch is achieved when the mass accretion rate 
approaches its Eddington limit (JGodon 1119971 ). the flow becomes spherically symmetric and 
the slim disk approximation for a one-dimensional BL becomes inaccurat e. Instead, the 



vertic ally averaged quantities represent the spherically symmetric quantities (JNarayan fc Yi 



19951 ). This implies that the results of the one-dimensional advective BL in the plane of the 
disk represent a massive spread layer engulflng the entire surface of th e star. Such simulation s 



have on l y been carried out for BL around FU Ori stars (e.g. in ID iPopham et al. I ( 1l993l ): 



Godon I (Il996bl ). and in 2D iKley fc Lin I (119961 )) and symbiotics (JGodon 



1996al ). 



For an accretion rate in the range 3 x 10^^^Mq/jt < M < 2 x IO^^Mq/jt, the matter 
does not spread far from the WD equator and most of the dissipated energy radiates back 
into the accretion disk. The BL remains confined to the plane of the disk, and one can 
expect the one-dimensional simulations to accurately represent the BL between the star and 
disk. In this case, the exa ct te mperature, siz e and rotation rate of the BL are given in 
Popham fc Naravan I (19951) an d I Go don et al. I ( 19951). The temperature is in the range ~ 
150. OOOK dCodon et al. Ill995h to a few 100,000K dPopham fc Naravan 1^9951 ). A modeling 
of the extreme UV (Chandra) spectra of s ome CVs (e.g. OY Car, SS Cyg fc WZ Sge) in this 
regime of M leads to such temperatures ( jMauche fc Raymond I l2000l : iMauche I |2004J ) . Here 
too, the BL emits in the soft X-ray regime with relatively more contribution to the extreme 
UV and far UV as the temperature decreases. The WD and the inner disk emit in the UV, 
the outer disk emits in the optical. The modeling we carry out in the present work is the 
modeling of such a BL. 

For M < 3 X 10"^°MQ/yr, the BL is optically thin and the only UV emission is coming 
from the inner disk edge irradiated by the optically thin BL. That inner edge forms a ring 
at a radius r ~ 1.2R^ — 1.6- R„,, rotating at Keplerian speed, with an effective temperature 
Teff ~ 100, OOOK-50, OOOK dNaravan fc Popham lEJ : IPopham Ill999h for an accretion rate 
of 3 X lO~^°M0/yr and 3 x lO~^^M0/yr, respectively. As the mass accretion rate decreases, the 
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ring radius increases and its effective temperature decreases. The BL emits in the hard X-ray 
range and the disk's inner edge ring emits in the UV if its temperature is high enough. Such 



UV e mission from the inner edge of the disk/ BL was identified in VW Hyi ( iPandel et al. 



20031 ) and in other dwarf novae in quiescence (jPandel et al. II2005I ). This UV contribution 



from the BL h as to be taken into a ccount when modehng the UV spectra of accreting WDs 
in quiescence ( JGodon fc Sion 1120051 ). The WD emits in the UV, as does the very inner disk. 
The rest of the disk emits in the optical. As the mass accretion rate decreases further, the 
disk eventually ceases to emit in the UV. 

As a consequence, the standard disk model is roughly correct for (say) r ^ 2i?^, and 
larger, at both high and low mass accretion rates. The exact value of r for which the 
standard disk model breaks down depends on the mass accretion rate and mass of the WD. 
At smaller radii the temperature in the disk-BL region increases rapidly compared to the 
standard model prediction. 



3. Data Analysis 



3.1. The FUSE Archival Data 



MV Lyr was observed with FUSE on May 6, 2003 (at 20h:07m:23s UT), just after it 
climbed into a high state that lasted for about 4yr. The Ft/S"^ exposure (D9050901) consists 
of 7637s of good exposure tira e. The data were processed with the latest and final version of 
Cal FUSE (IDixon et al. II2007I) following the same procedure as in previous works (for details 



see 



Godon et al. 



(I2OO9I)). 



The FC/^E^ spectrum of MV Lyr (Figure 1) is characterized by broad and deep absorption 
lines from highly ionized species and the absence of Hydrogen Lyman lines (except for sharp 
lines from the ISM). These absorption lines are Niv (~923A ), Svi (933.5A & 944.5A ), 
Ovi (1131.6 & II37.3A ), Siv (1062.6 & 1073A ), Sim (~ llOSA ), Pv (lllsA ), Siiv 
(1122.5 & 1128. 3A ), and Cm (1175A ). Such deep and broad absorption lines are often 
observed in CVs in high state, even in systems with a large inclination where the Keplerian 
velocity broadening is expe cted to produce a ra ther smooth continuum (e.g. see the FUSE 
spectrum of V3885 Sgr in iLinnell et al. I (120091 )). These broad and deep absorption lines 
most likely form in a hot corona above the heated stellar surface, BL or disk. Additional 
transition wavelengths have been marked on the spectrum in Figure 1 and can be associated 
with some spectral features; however, they are not clearly detected. 

There are a few sharp absorption lines from the interstellar medium such as Si ll 1020A, 
Cll 1036A, An 1048.2 & 1066. 7A, and Fell 1145A. The sharp emission lines are artifacts 
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due to air glow or direct reflection of sunlight inside the telescope (e.g., such as the Cm 
977A & Hell 1168A lines). 



3.2. TLUSTY and SYNSPEC Codes 

The modeling of the FUSE spectiuTCi is carried out in steps. First, spherically symmetric 
stellar atmosphere struct ure as well as the vertical structure of disk rings are computed 



using the code TLUSTY (iHubeny Ill988l ). For the stellar structure the input parameters are 
the surface gravity, effective temperature and the surface composition of the star. For the 
composition it is often enough to take Hydrogen and Helium and neglect the metals as long 
their abundance is low (say similar to solar composition), in most cases this does not affect 
the computed atmospheric structure. For high temperatures we use the NLTE option. For 
the disk rings the input is the local mass accretion rate, mass of the accreting star, radius of 
the star, and radius of the ring and solar composition. Alternatively, one can input for each 
ring the effective temperature, effective vertical gravity and the column mass at midplane. 

Second, the code SYNSPEC is used to derive the detailed radiation and flux distribution 



of continuum and lines (JHubeny fc Lanz Ill995l ). SYNSPEC generates the output spectrum 



of the stellar atmosphere structure and disk vertical structure computed by TLUSTY. The 
input for SYNSPEC is the same as for TLUSTY and it uses the structure computed by 
TLUSTY as an input. Composition is specified here with an additional input file, and can 
include species up to iron or even zinc. In the present modeling we included only H, He, C, 
N, O, Si, P, and S. 

In the next step the code ROTIN is used to account for rotational and instrumen- 
tal broadening of the l ines for the stellar synthetic spectrum, while the code DISKSYN 



(IWade &: Hubeny Ill998l ) is used to combine the disk rings together (where we substitute the 
inner rings with hot BL rings if a BL is computed) and account for rotational broadening 
due to Keplerian motion. Limb darkening is also included in the codes DISKSYN & ROTIN. 

Boundary layer synthetic spectra are generated by computing inner disk rings with an 
elevated temperature matching the BL temperature obtained from theoretical estimates. 
The rings of the BL are then substituted to the inner disk rings. We then use a reduced xl 
fitting technique to find the best fit to the observed spectrum. During the xt fitting process 
of the FUSE spectrum taken in the high state the mass of the WD, the inclination of the 
system, and the composition are kept constant because they have been firmly established 
from previous observations. The remaining parameters, i.e. the temperature of the WD, the 
mass accretion rate, the number of BL rings included, and the temperature of the BL rings 
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are allowed to vary as they are unknown. The distance to the system is obtained a posteriori 
from the spectral fitting as an output parameter; resulting models with a distance that does 
not agreed with the assumed distance are rejected (see next section). 

In the present work we use the following versions of the software: TLUSTY202, SYN- 
SPEC48, & R0TIN4 for WD photospheric spectra runing on a Linux platform (Cygwin-X) 
with a GNU FORTRAN 77 compiler; TLDISK195 (a previous variant of TLUSTY), SYN- 
SPEC43, & R0TIN3 and DISKSYN7 for disk and BL spectra runing on a Unix machine also 
using a GNU FORTRAN 77 compiler (the GNU FORTRAN 77 compiler runs slightly dif- 
ferently under Unix and Linux operating systems). The software packages and user's guides 
are available online (free download) at ,http://nova.astro. umd.edu/, . 



Results 



In the following spectral modeling of the FUSE spectrum of MV Lyr in its high state, as 
we derive the mass accretion rate and characterize the BL of the system, we wish to use the 
mass of the WD and the distance to the system as given. These parame t ers w ere derived from 
the FUSE spectrum of the system in the low state by iHoard et al. I (120041 ). who obtained 
a WD surface temperature T = 47, OOOK, a gravity logs' = 8.25, non-solar composition 
Z ^ 0.3^0, and a distance d ~ 505ifc 50pc. Since these par ameters are extremely important, 
we decided to confirm the results of iHoard et al. I (120041 ) and we carried out a modeling 
of the FUSE spectrum of MV Lyr obtained in the low state. We obtained a WD surface 
temperature T = 44, OOOK, a gravity \ogg = 8.1, a projected stellar rotational velocity 
Vrot^mi = 200km/s, non-solar composition Z ^ 0.5Zq, and a distance d ~ 550 ± 50pc. The 
small difference in the parameters derived is du e to the facts that ( i) iHoard et al. I (J2004J ) 



did not in clude rotational broadening ; and (ii) iHoard et al. I ( 120041 ) used the mass-radius 
relation of iHamada fc Salpeter I ( Il96ll ) for zero-temperature WDs ( giving M^a = 0.73Mq 
for logs' = 8.25), while we used the C/O-core WD model of IWood I (Il995l ) for a ~45,000K 
WD. (giving M^a = 0.73Mq for logs = 8.1 but M^a ~ O.8M0 for logs = 8.25), We 
note, however, that our resu l ts are in complete agreement with the modeling carried out by 
E.M.S. within iHoard et al. I (120041 ) 's work. Therefore, in the following we use the following 
parameters: M^d = O.8M0, d ~ 500 — 550pc, T^d = 44, OOOK and above. The abundances 
and stellar rotational velocity have no effects on the fitting to the FUSE spectrum in the high 
state as it is dominated by emission from the solar composition disk with a large Keplerian 
velocity broadening. 

In the fitting of the high state spectrum, no attempt is made to model the broad and 
deep absorption lines. These lines form in a hot corona above the disk/BL and cannot be 
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modeled with the synthetic spectral code TLUSTY/SYNSPEC. The code models the lines 
created in the disk itself, but not in a hot corona above it. As a consequence, only the 
observed continuum is considered and the broad absorption lines are masked. 

The WD model used initially in the combined BL + disk + WD fits to the FUSE spec- 
trum of MV Lyr in the high state has T = 44, OOOK. In the high state the WD temperature 
could be elevated due to continuing accretion. Therefore, models with larger WD tempera- 
ture are also computed (50, OOOK and up) but they give very similar results. This is due to 
the fact that the WD is not the main emitting component of the UV spectrum. The details 



are given below. We assume here i = 10° ± 3° (ISchneider et al. Ill98ll : ISkillman et al. Ill995 



Linnell et al. Il2005[ ) 



The first four models presented in Table 2 are the 44, OOOK WD -|- standard disk models 
with increasing mass accretion rate. All these models do not have enough flux in the shorter 
wavelength range to fit the FUSE spectrum. The model with M = 5 x lO~^M0/yr gives a 
distance of 666pc, too large to be acceptable, and with M = 1 x 10^^Mq/jt the distance 
becomes twice the value found in the low state. As an example, a model with M = 3.5 x 
lO^^M0/yr is shown in Figure 1, giving a distance of 557pc. The flux deficiency in the shorter 
wavelengths is clearly seen. The inclusion of a hotter WD does not provide a significant 
improvement of the model fit. 

In order to increase the flux in the shorter wavelengths, the temperatures of the two 
inner rings of the standard disk model are modified to represent the BL . The first ring is 
located at ri = 1.05i?* and the second is at r2 = 1.20-R*. The temperatures of the rings (Ti 
& T2 respectively) are listed in Table 2. In the standard disk model these temperatures are 
below 50, OOOK for the accretion rate considered here. For the modeling of the BL , rings are 
compute d with temperatures between 100, OOOK and 175, OOOK, in agreemen t with the BL 
models of Icodon et al. I f llQQsl ) (~ 125, OOOK) and lPopham fc Naravan Nl995l ) (~ 180, OOOK) 
for the mass accretion rates considered here. In doing so, one is able to construct two models 
of BL: thin (first ring only) and extended (first two rings). 

The standard disk model with a mass accretion rate of 3 x 10^^Mq/jt fits the long 
wavelength end of the FUSE spectrum using the assumed distance to MV Lyr. For this 
reason models with M > 3 x lO~^M0/yr are not considered when including the BL, since 
the BL increases the flux of the model and therefore its distance becomes too large. Also, 
since MV Lyr is in a high state, mass accretion rates below 1 x 10^^Mq/jt are not considered 
(these points are discussed in the section 5). 

First, an accretion disk with a mass accretion rate of 1 x lO~^M0/yr is considered, and 
the BL temperature and size are then varied. It is found that the best model flts lead to a 
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distance of only 350pc-400pc. Since the mass accretion rate is fixed to 1 x IO^^Mq/jt, the 
only way to increase the distance is to increase the contribution of the BL , by increasing its 
temperature and/or size. In doing so, the models that give an acceptable distance (of say at 
least ~ 430pc) have actually too much flux in the shorter wavelengths, a sign that the BL 
contributes too much flux. Only a few of these models are listed in Table 2. These models 
are not better than the standard disk models: their xt i^ ^s large and/or their distance is 
too short. The best fit models are obtained for a thin (one ring) BL with a temperature 
of ~150,000K-175,000K. The distance for these models is, however, far too short. The 
inclusion of a heated WD does not improve the models and produces only a small increase 
in the distance. 

Next, to obtain models with a larger distance, the mass accretion rate is increased to 
2 X 10~^MQ/yr. This has the effect of decreasing the relative flux contributed by the BL. 
These models agree with the assumed distance and have a lower x^, especially the two-ring 
BL models. The best one-ring model is for a BL temperature of 175,000K. This model has 
a relatively low xt ^^^ agrees well with the distance. For the two-ring BL models, there is 
little difference (in terms of xt) between a temperature of 100,000K, 125,000K and 150,000K 
(or any combination of those) as they produce x^ ~ 1 with d ~ 500 ± 50pc. A 150,000K 
extended (2 rings) BL model with disk and WD is shown in Figure 2. This model fits the 
flux in the shorter wavelength range very well (A < 970A , ignoring the broad absorptions 
of N IV & S Vl) but produces too much flux around 980-lOlOA . Extended BL models with 
a lower temperature (e.g. 100,000K) do not fit the shorter wavelengths as well, but on the 
other side they better fit the 970-lOloA region and, consequently, they have the same xl- 
In a last effort to improve the modeling, the temperature of the WD is increased to see how 
it affects the results. As the WD temperature increases the hotter BL models deteriorate 
slightly while the cooler BL models improve slightly. These models are listed in the lower 
part of Table 2 A fit with a 50,000K WD, an extended 100,000K BL and a 2 x lO-^Mo/yr 
accretion disk model is shown in Figure 3. 



5. Discussion and Conclusions 

Overall, the best fit to the FUSE spectrum of MV Lyr in its high state consists of a 
composite WD + disk -|- BL model with a mass accretion rate M = 2 x 10~^MQ/yr. There 
are two equally acceptable BL solutions, either a broad extended (two-ring) BL with a rather 
low temperature (in the range 100,000K-150,000K), or a thin (one-ring) BL with a higher 
temperature (175,000K). The WD temperature is probably higher than the 44,000K found 
in the low state, due to ongoing accretion but it affects only slightly the quality of the fitting 
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of the models, we discuss this further at the end of this section. 

Since the modehng of the FUSE spectrum generates more than one solution over a 
finite spectral range, one must consider the bolometric luminosity of each component for 
comparison with equations (2) & (6). For each model, the total luminosity of the BL is 
computed assuming it radiates as a black body. This is a reasonable assumption, as in this 



regime the BL is most probably optically thick (JGodon et al. I Il995t iPopham fc Narayan 



19951 ). The two components of temperature Ti & T2 (given in Table 2), and with radii 
Ti = 1.05R^ and r2 = 1.2i?* are added. One then simply uses the Stefan-Boltzmann law and 
area of each ring to find the total BL luminosity. The bolometric luminosity obtained is then 
compared with the theoretical expectation in eq.(6) using eq.(2). For almost all the disk+BL 
models listed in Table 2 it is found that the computed bolometric BL luminosity is larger than 
the disk luminosity. The discrepancy is even larger for the model with M = 1 x lO^^M0/yr. 
The only models for which the BL bolometric luminosity is smaller than the disk luminosity 
(eq.6) is for an extended BL (two-ring model) with a temperature as low as 100,000K. For 
that model the bolometric BL luminosity is 93% of the disk luminosity, implying a rotational 
velocity (eq.6) of less than 4% its Keplerian value or about 135km/s. 

With an inclination i = 10° ± 3° and a projected rotational velocity of 200 ± 50km/s 
(as derived in the low state), the (non-projected) rotational velocity should be in the range 
667km/s< Vrot <2050km/s, or 0.18 < Q^/Qk < 0.54 (the lower limit is for a velocity of 
150km/s and i = 13°, while the upper limit is for 250km/s and i = 7°). This produces a BL 
luminosity in the range 0.68 > Lbl/ Ldisk > 0.21. For the BL model to agree with the upper 
limit, one could reduce the BL temperature (by 10,000K) or increase the mass accretion rate 
(by 1/3, such a model is listed at the end of Table 2) or both. However, one notes that the 
error in the WD mass of O.IMq (~ 0.2 in logg) produces a computed Lbl/ Ldisk between 
0.55 (for a O.9M0 WD mass) and ~1.5 (for a 0.7Mq WD mass). Since the BL model used 
here is rather simplistic, a perfect agreement between the computed bolometric luminosity 
of the BL and that computed from eq.6 is not expected. However, within the limits of the 
errors, the results of the BL model are consistent with a broad BL (of size ~ 0.2i?*) with a 
rather low temperature (~ 100, OOOK or slightly less) and a mass accretion rate of the order 
of ~ 2 X lO~^M0/yr (or slightly larger). These models are located the lower part of Table 2 
and the fit improves (lower xt ^''^d better distance) as the temperature of the WD increases, 
reaching a best fit for T = 70, OOOK. Though the temperature of the WD is certainly not 
that high, this certainly points to the fact that the temperature of the WD is elevated during 
the high state reaching ~ 50, OOOK or a little higher. 
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Table 1. System Parameters for MV Lyr 



Parameter 



Value 



Reference 



(1995} 



this work; Hoard et al. (2004) 
based on models of Wood (1995) 
Hoard et al. ^2004) 
Schneider et al. (1981J; Skillman et al. 
Linnell et al. ('20051 
Bruch fc Engel ('19941 
Skillman et al. ('19951 
this work; Hoard et al. (2004) 
2 X lO-^Me/yr This work; Linnell et al. (2005) 



M^d 


O.73-O.8M0 


J^wd 


7,440km 


M2nd 


0.3Mq 


i 


10° ±3° 


E(B-V) 





P 


3.19hr 


d 


~ 500 - 550p 


M 


2 X lO-^MfT,/- 
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Table 2. Accretion Disk & Boundary Layer Synthetic Spectral Model Fits 



Model 


Mx 10^ 


Twd 


d 


Ti 


T2 


WD/disk 


xl Figure 




(Mo/yr) 


(lO^K) 


{pc) 


(lO^K) 


(lO^k) 


% 




disk 


1.0 


44 


301 


25.9 


32.2 


19/81 


3.066 


disk 


2.0 


44 


428 


30.7 


38.3 


10/90 


1.923 


disk 


3.5 


44 


557 


35.4 


44.1 


6/94 


1.379 1 


disk 


5.0 


44 


666 


38.7 


48.2 


4/96 


1.270 


disk+BL d 


1.0 


44 


370 


175 


32.2 


13/87 


0.856 


disk+BL 1 


1.0 


44 


349 


150 


32.2 


14/86 


1.083 


disk+BL 2 


1.0 


44 


336 


125 


32.2 


16/84 


1.376 


disk+BL 3 


1.0 


44 


327 


100 


32.2 


16/84 


1.376 


disk+BL 4 


1.0 


44 


448 


150 


150 


9/91 


1.302 


disk+BL d 


2.0 


44 


480 


175 


38.3 


8/92 


1.019 


disk+BL 1 


2.0 


44 


464 


150 


38.3 


8/92 


1.203 


disk+BL 2 


2.0 


44 


454 


125 


38.3 


9/91 


1.358 


disk+BL 3 


2.0 


44 


447 


100 


38.3 


9/91 


1.485 


disk+BL 4 


2.0 


44 


538 


150 


150 


6/94 


1.058 2 


disk+BL 5 


2.0 


44 


507 


125 


125 


7/93 


0.993 


disk+BL 6 


2.0 


44 


484 


100 


100 


8/92 


1.056 


disk+BL da 


2.0 


44 


516 


175 


100 


7/93 


0.985 


disk+BL 7 


2.0 


44 


500 


150 


100 


7/93 


1.000 


disk+BL 8 


2.0 


44 


491 


125 


100 


7/93 


1.025 


disk+BL da 


2.0 


44 


516 


175 


100 


7/93 


0.985 


disk+BL da 


2.0 


50 


521 


175 


100 


8/92 


0.995 


disk+BL 4 


2.0 


44 


538 


150 


150 


6/94 


1.058 2 


disk+BL 4 


2.0 


50 


542 


150 


150 


8/92 


1.078 


disk+BL 4 


2.0 


60 


548 


150 


150 


10/90 


1.082 


disk+BL 4 


2.0 


70 


555 


150 


150 


12/88 


1.090 


disk+BL 6 


2.0 


44 


484 


100 


100 


8/92 


1.056 


disk+BL 6 


2.0 


50 


489 


100 


100 


9/91 


1.030 3 
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Table 2 — Continued 



Model 


Mx 109 

{Mrjyr) 


Twd 
(lO^K) 


d 

(pc) 


Ti 
(lO^K) 


T2 

(lO^k) 


WD/disk 

% 


xl 


Figure 


disk+BL 6 


2.0 


60 


496 


100 


100 


12/88 


1.012 




disk+BL 6 


2.0 


70 


503 


100 


100 


14/86 


0.989 




disk+BL 6 


2.7 


50 


545 


100 


100 


8/92 


1.068 
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Fig. 1. — The FUSE spectrum of MV Lyr in its high state is modeled with a 44,000K 
WD (dotted hne), with a mass 0.8Mq and a standard accretion disk (dashed hne) with 
M = 3.5 X lO~^M0/yr and i = 10°. The combined model WD + disk is the solid black 
line. This gives a matching distance of 557pc. The model cannot match the flux level in the 
short wavelength range (< 950A). In order to fit that region the mass accretion rate has to 
be increased, and this increases the distance to twice the distance of MV Lyr. 
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Fig. 2. — The FUSE spectrum of MV Lyr in its high state is modeled with a composite 
WD + disk model which includes the BL . The WD model (dotted line) has a temperature 
of 44,000K and a mass O.SMq, the accretion disk model (dashed line) has a mass accretion 
rate of M = 2.0 x 10~^Mq/jt and i = 10°. The BL is modeled as the two inner rings of 
the disk with a temperature of 150,000K. The combined model WD + disk is the solid black 
line. This gives a distance of 538pc. The WD contributes only 6% of the flux in the FUSE 
spectral range. 
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Fig. 3. — The FUSE spectrum of MV Lyr in its high state is modeled with a composite WD 
+ disk model which includes the BL . The WD model (dotted line) has a temperature of 
50,000K and a mass O.SM©, the accretion disk model (dashed line) has a mass accretion rate 
of M = 2.0 X IO'^Mq/jt and i = 10°. The BL is modeled as the two inner rings of the disk 
with a temperature of 100,000K. The combined model WD + disk is the solid black line. 
This gives a distance of 489pc. The WD contributes 9% of the flux in the FUSE spectral 
range. 



